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ABSTRACT 

We  present  results  from  an  ongoing  multiwavelength  radial  velocity  (RV)  survey  of  the  Taurus-Auriga  star-forming 
region  as  part  of  our  effort  to  identify  pre-main-sequence  giant  planet  hosts.  These  1-3  Myr  old  T  Tauri  stars  present 
significant  challenges  to  traditional  RV  surveys.  The  presence  of  strong  magnetic  fields  gives  rise  to  large,  cool 
star  spots.  These  spots  introduce  significant  RV  jitter  which  can  mimic  the  velocity  modulation  from  a  planet-mass 
companion.  To  distinguish  between  spot-induced  and  planet-induced  RV  modulation,  we  conduct  observations  at 
~6700  A  and  ~2.3  //in  and  measure  the  wavelength  dependence  (if  any)  in  the  RV  amplitude.  CSHELL  observations 
of  the  known  exoplanet  host  G1  86  demonstrate  our  ability  to  detect  not  only  hot  Jupiters  in  the  near-infrared  but 
also  secular  trends  from  more  distant  companions.  Observations  of  nine  very  young  stars  reveal  a  typical  reduction 
in  RV  amplitude  at  the  longer  wavelengths  by  a  factor  of  ~2-3.  While  we  cannot  confirm  the  presence  of  planets 
in  this  sample,  three  targets  show  different  periodicities  in  the  two  wavelength  regions.  This  suggests  different 
physical  mechanisms  underlying  the  optical  and  the  V-band  variability. 
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1.  INTRODUCTION 

The  discovery  of  over  760  exoplanets8  in  the  past  twenty  years 
has  revealed  that  planetary  systems  are  common  and  diverse. 
Pulsar  planets  (Wolszczan  1994),  hot  Jupiters  (Mayor  &  Queloz 
1995),  and  retrograde  orbits  (Hebrard  et  al.  2011)  repeatedly 
challenge  our  assumptions  about  planet  formation.  After  two 
decades  of  exoplanet  discoveries,  the  processes  underlying 
planet  formation  remain  unclear.  Lacking  direct  observational 
inputs,  theorists  must  deduce  formation  mechanisms  from 
observations  of  mature  systems.  A  more  preferable  approach  is 
to  catalog  the  planet  population  around  very  young  (~  1  Myr  old) 
stars  while  these  nascent  worlds  still  reside  in  their  formation 
environments. 

Unfortunately,  observations  of  pre-main-sequence  (PMS) 
stars  are  complicated  by  star  spots,  jets,  accretion,  and  cir- 
cumstellar  disks.  The  most  successful  attempts  at  young  planet 
detection  to  date  have  come  from  direct  imaging  surveys.  For 
example,  Lafreniere  et  al.  (2010)  presented  evidence  for  a  gi¬ 
ant  planet  around  the  ~5  Myr  old  Upper  Sco  solar  analogue 
1RXS  J 160929. 1-210524.  They  report  the  presence  of  an 
~8  Mjup  companion  at  a  distance  of  330  AU  based  on  JHK 
adaptive  optics  images  from  NIRI  on  Gemini  North.  Another 
direct  imaging  detection  is  reported  in  Kraus  &  Ireland  (2012). 
They  present  evidence  for  a  faint  point  source  ~16  AU  from 
the  2  Myr  old  LkCa  1 5  using  non-redundant  aperture  masking 
interferometry  on  Keck  II.  They  interpret  this  source  as  a  ~6 
Mjup  protoplanet.  These  observations,  and  others  like  them,  are 
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providing  the  first  glimpse  into  the  mechanisms  of  giant  planet 
formation. 

There  are,  however,  limitations  to  what  can  currently  be 
accomplished  with  direct  imaging.  This  technique  can  only 
detect  massive  planets  at  wide  separation  from  their  host  stars. 
Furthermore,  imaging  does  not  provide  a  direct  measure  of 
companion  mass;  the  derived  masses  are  model  dependent. 
Radial  velocity  (RV)  surveys,  which  have  confirmed  ~95%  of 
the  known  exoplanet  population,9  obtain  mass  limits  and  are 
more  suited  to  identifying  planets  on  close  orbits.  In  particular, 
“hot  Jupiters”  are  readily  accessible  via  RV  observations. 

Low-mass  young  stars,  however,  are  challenging  targets  for 
traditional  RV  surveys.  These  targets  are  faint  at  optical  wave¬ 
lengths  owing  to  late  spectral  types,  large  distances  (>100  pc), 
and  extinction  from  natal  dust  clouds.  Young  stars  also  have 
strong  magnetic  fields  (e.g.,  Johns-Krull  2007)  that  generate 
large,  cool  star  spots.  Spots  impact  RV  surveys  of  young  stars 
by  introducing  significant  jitter  which  can  mimic  the  RV  mod¬ 
ulation  induced  by  a  planet  (Saar  &  Donahue  1997;  Hatzes 
2002;  Desort  et  al.  2007;  Reiners  et  al.  2010).  Recent  attempts 
at  detecting  substellar  companions  in  young  stellar  populations 
(10-100  Myr)  via  RV  monitoring  have  generally  been  unsuc¬ 
cessful  (e.g.,  Paulson  et  al.  2004;  Paulson  &  Yelda  2006),  likely 
because  of  small  sample  sizes  and  the  intrinsic  RV  variability 
of  the  targets.  For  very  young,  1-3  Myr  old  stars,  the  RV  noise 
problem  is  even  more  acute.  Recently,  van  Eyken  et  al.  (2012)  re¬ 
ported  observations  of  a  possible  hot  Jupiter  orbiting  a  2-3  Myr 
old  weak-lined  T  Tauri  star  in  Orion  ( P  0.45  days).  This 
interpretation  is  based  on  transit  data  from  the  Palomar  Tran¬ 
sit  Factory  project  in  conjunction  with  optical  RV  observations 
from  the  Hobby-Eberly  Telescope  and  Keck  I.  While  this  is  an 
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intriguing  result,  their  data  are  dominated  by  spot  modulation. 
Confirmation  must  wait  for  further  observations,  particularly  in 
the  near-infrared  (N1R).  The  youngest  confirmed  RV  planet  to 
date  is  a  ~6  Mjup  planet  on  an  850  day  orbit  around  the  100  Myr 
old  star  HD  70573  (Setiawan  et  al.  2007). 

The  success  of  a  young  star  RV  survey  depends  on  its 
ability  to  distinguish  between  spot-induced  and  companion- 
induced  RV  modulation.  Spot-induced  RV  variations  in  young 
stars  have  been  identified  by  correlations  between  RVs  and 
spectral  line  bisector  spans  (e.g.,  Hatzes  et  al.  1997);  planet 
hosts  do  not  exhibit  any  such  correlation.  However,  there 
may  also  be  no  correlation  if  the  absorption  lines  are  not 
spectroscopically  resolved,  i.e.,  the  projected  rotation  velocity 
( v  sin  i)  of  the  star  is  comparable  to  or  less  than  the  velocity 
resolution  of  the  spectrograph  (Desort  et  al.  2007;  Huelamo 
et  al.  2008;  Prato  et  al.  2008).  A  potentially  more  reliable 
method  for  distinguishing  between  spots  and  planets  leverages 
the  wavelength  dependence  of  a  spot-induced  RV  modulation 
amplitude. 

The  reflex  motion  caused  by  a  planet  affects  all  wavelengths 
equally.  However,  the  amplitude  of  any  spot-induced  RV  vari¬ 
ability  will  be  smaller  at  longer  wavelengths  (e.g.,  Huelamo 
et  al.  2008;  Prato  et  al.  2008;  Reiners  et  al.  2010;  Ma  &  Ge 
2012).  This  wavelength  dependence  arises  because  of  the  flux- 
temperature  scaling  in  the  Rayleigh-Jeans  limit  of  blackbody 
radiation;  the  contrast  between  a  photosphere  and  a  cooler  star 
spot  decreases  at  longer  wavelengths  (e.g.,  Vrba  et  al.  1986; 
Carpenter  et  al.  2001).  Therefore,  observations  in  the  visible 
and  NIR  may  distinguish  between  stellar  activity  and  a  true 
companion  by  comparing  the  RV  amplitudes  at  the  two  wave¬ 
lengths.  Furthermore,  the  late  spectral  types  of  T  Tauri  stars  pro¬ 
duce  SEDs  with  peak  emission  around  1-2  /xm  and  extinction 
of  these  partly  embedded  sources  is  lessened  at  longer  wave¬ 
lengths.  NIR  observations  can  therefore  lead  to  an  increased 
signal-to-noise  ratio  (S/N)  and  improved  RV  precision. 

NIR  spectroscopy  has  only  recently  begun  to  play  a  role 
in  the  search  for  substellar  companions.  Martin  et  al.  (2006) 
combined  optical  and  //-band  observations  to  look  for  giant 
planets  around  the  young  brown  dwarf  LP  944-20.  They 
concluded  that  the  observed  optical  RV  modulations  were  driven 
by  inhomogeneous  surface  features  (i.e.,  clouds).  Blake  et  al. 
(2007)  used  high-resolution  K- band  spectroscopy  to  investigate 
the  presence  of  giant  planets  around  a  population  of  L  dwarfs, 
achieving  a  precision  of  300  m  s'1 .  They  found  no  evidence 
for  companions  with  M  sin  i  >  2  Mjup  and  P  <  3  days  in  their 
sample  of  nine  targets.  Setiawan  et  al.  (2008)  reported  the  optical 
detection  of  a  ~  10  Mjup  planet  on  a  3.5  day  orbit  around  the 
10  Myr  old  star  TW  Hydra.  However,  //-band  observations  by 
Huelamo  et  al.  (2008)  revealed  a  strong  wavelength  dependence 
in  the  velocity  amplitude.  The  low  H- band  variability  casts  doubt 
on  the  presence  of  a  companion  and  suggests  that  spots  cause 
TW  Hydra’s  optical  RV  variations.  Prato  et  al.  (2008)  observed 
DN  Tau  and  V836  Tau — two  potential  planet-host  candidates 
selected  on  the  basis  of  their  optical  variability — in  the  K  band. 
To  within  their  measurement  precision,  they  detected  no  NIR 
RV  modulation,  implying  that  spots  produce  the  apparent  optical 
RV  variability. 

Interest  in  targeting  M  dwarfs  as  hosts  of  habitable  planets 
has  spurred  vast  improvements  in  NIR  precision,  bringing  it 
closer  to  that  of  optical  surveys.  Bean  et  al.  (2010)  report  a 
long-term  precision  in  the  H  band  of  ~5  m  s~  1  on  late-M 
dwarfs  using  an  ammonia  gas  cell  on  CRIRES  at  the  Very  Large 
Telescope  (VLT).  Figueira  et  al.  (2010)  report  a  comparable 
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precision,  also  on  CRIRES,  but  using  telluric  absorption  features 
as  a  wavelength  reference.  Based  on  six  years  of  NIRSPEC 
data,  Blake  et  al.  (2010)  report  a  precision  of  50  m  s~  1  for 
their  sample  of  K  dwarfs  and  200  m  s-1  for  L  dwarfs,  also 
through  the  use  of  telluric  lines.  Bailey  et  al.  (2012)  also 
report  50  m  s  1  precision  on  mid-M  dwarfs,  increasing  to 
~80-170  m  s-1  for  young  (~  10  Myr  old)  stars,  using  NIRSPEC 
with  telluric-based  wavelength  calibration.  While  these  results 
are  encouraging,  these  efforts  are  all  focused  on  large  (8-10  m) 
aperture  telescopes.  Very  little  work  has  been  done  to  develop 
precision  NIR  RV  techniques  for  smaller  aperture  telescopes 
where  more  observing  time  is  available  to  the  community.  We 
recently  presented  our  technique  for  achieving  ~50  m  s  1  on 
late-K  and  early-M  dwarfs  using  the  CSHELL  NIR  spectrograph 
on  the  Infrared  Telescope  Facility  (IRTF;  Crockett  et  al.  2011). 
Anglada-Escude  et  al.  (2012)  also  show  very  promising  results 
with  their  successful  efforts  at  obtaining  20-30  m  s”1  precision 
in  the  K  band,  also  on  CSHELL,  using  a  methane  isotopologue 
gas  cell. 

Since  2004,  we  have  been  conducting  an  RV  survey  in  the 
Taurus-Auriga  star-forming  region  in  an  attempt  to  catalog  the 
giant  planet  and  brown  dwarf  population  around  very  young 
stars  (Huerta  et  al.  2008;  Prato  et  al.  2008;  Crockett  et  al.  2011; 
Mahmud  et  al.  2011).  In  this  paper,  we  present  new  and  updated 
analyses  on  nine  targets  including  known  spotted  stars  and 
potential  young  planet  hosts;  the  entire  sample  will  be  discussed 
in  a  later  paper  once  the  entire  survey  is  complete.  In  Section  2, 
we  provide  an  overview  of  our  observing  strategy.  Our  data 
reduction  pipeline  and  RV  analysis  techniques  are  presented  in 
Section  3,  including  results  from  RV  standards  used  to  quantify 
our  long-term  precision.  In  Section  4,  we  discuss  the  results 
of  these  analyses  for  each  of  our  T  Tauri  targets,  as  well  as 
updated  results  on  the  known  exoplanet  host  G1  86.  Finally,  in 
Section  5,  we  discuss  the  relevance  of  these  results  and  how 
they  can  inform  future  studies  of  the  young  planet  population. 

2.  OBSERVATIONS 

Our  target  list  mainly  consists  of  very  young  (1-3  Myr) 
T  Tauri  stars  in  the  nearby  (~140  pc;  Kenyon  et  al.  2008) 
Taurus-Auriga  star-forming  region,  many  of  which  we  have 
been  observing  for  several  years.  These  targets  are  selected 
from  Herbig  &  Bell  (1988)  based  on  the  following  criteria:  (1) 
V  <  15,  (2)  v  sin  i  <  20  km  s"  1  to  ensure  relatively  narrow,  deep 
photospheric  lines  for  accurate  RV  analysis,  and  (3)  GKM  stars 
to  maximize  S /N  in  the  red  part  of  the  spectrum  for  our  infrared 
observations.  PMS  stars  of  late-K  and  early-M  spectral  types 
will  evolve  to  become  late-G,  early-K  solar  analogues.  They  of¬ 
fer  a  good  compromise  between  having  the  deep  photospheric 
CO  lines  needed  for  NIR  RV  measurements  (Section  3.3)  and 
being  massive  enough  to  offer  a  reasonable  chance  of  hosting 
a  giant  planet  (Johnson  et  al.  2010).  Of  the  129  T  Tauri  stars, 
~35%  are  classical  (actively  accreting)  and  ~65%  are  weak- 
lined  (little  to  no  accretion).  We  also  include  14  targets  from  the 
older  (~  100  Myr)  Pleiades  open  cluster.  About  20%  of  our  sam¬ 
ple  are  known  to  be  close  (>  O'.' 05)  binaries  which  are  unresolved 
in  our  observations.  Although  these  targets  are  RV  variables,  the 
timescale  is  on  the  order  of  tens  to  hundreds  of  years  and  thus 
does  not  impact  our  search  for  short-period  variability. 

We  conduct  observations  in  both  the  optical  (~6200-7100  A) 
and  the  K  band  (~2.3  /xm)  to  identify  any  wavelength  depen¬ 
dence  on  RV  variability.  Our  observing  strategy  starts  with  ob¬ 
taining  optical  observations  of  a  subset  of  targets  every  night 
for  roughly  a  week  at  a  time.  A  week-long  observing  window 
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matches  the  typical  rotation  period  of  T  Tauri  stars  and  falls 
within  the  range  of  companion  orbital  periods  that  we  are  capa¬ 
ble  of  detecting.  On  the  basis  of  the  optical  data,  we  then  group 
our  observed  targets  into  two  categories:  (1)  stars  that  exhibit 
correlations  between  RV  variations  and  line  bisector  changes 
and  (2)  stars  that  show  no  such  correlations.  The  RV  variations 
of  stars  in  Group  1  are  likely  attributable  to  star  spots  (Huerta 
et  al.  2008),  while  stars  in  Group  2  may  harbor  substellar  com¬ 
panions.  However,  since  a  lack  of  RV-bisector  correlation  is 
not  sufficient  for  definite  identification  of  a  companion,  we  test 
the  planet  hypothesis  by  conducting  follow-up  observations  of 
Group  2  stars  in  the  NIR  (Prato  et  al.  2008;  Mahmud  et  al.  201 1). 

With  this  strategy,  we  only  target  stars  that  show  no  correlation 
between  RV  and  bisector  span  for  further  observations  in  the 
NIR.  The  question  then  arises  as  to  how  many  RV-bisector 
pairs  are  needed  to  identify  a  correlation  between  these  two 
quantities  for  a  typical  spotted  T  Tauri  star.  To  address  this 
question,  we  examined  randomly  sampled  subsets  of  various 
sizes  from  the  RV  and  bisector  data  of  the  spotted  young  star 
LkCa  19  (Huerta  et  al.  2008).  Using  Monte  Carlo  methods,  we 
found  a  correlation  between  the  RV  and  bisector  data  90%  of  the 
time  at  a  false  alarm  probability  (FAP)  of  0.05  when  we  sampled 
12  observations.  Hence,  with  12  measurements  per  star,  1  in  20 
legitimate  exoplanet  systems  will  be  erroneously  rejected  by 
this  test  and  about  1  in  10  spotted  stars  with  clear  correlations 
will  mistakenly  pass  through  to  the  next  stage  of  observations. 
While  more  observations  per  star  would  improve  these  numbers, 
12  observations  strike  a  reasonable  balance  between  the  need  to 
have  a  large  sample  of  stars  and  good  time-coverage  per  star. 

2.1.  McDonald  Observatory 

Visible  light  spectra  were  taken  at  the  McDonald  Observatory 
2.7  m  Harlan  J.  Smith  telescope  on  the  Robert  G.  Tull  Coude 
Spectrograph  (Tull  et  al.  1995).  The  cross-dispersed  coude 
covers  a  wavelength  range  of  3986-9952  A.  We  isolated  the 
orders  with  strongest  S/N,  6233-7109  A,  for  RV  analysis.  A 
172  slit  yielded  a  spectral  resolving  power  of  R  =  (7, /A/,)  ~ 
60,000.  Integration  times  were  typically  1800  s  (depending  on 
conditions)  and  typical  seeing  was  ~2".  We  took  thorium-argon 
(ThAr)  lamp  exposures  before  and  after  each  target  observation 
for  wavelength  calibration;  typical  rms  values  for  the  dispersion 
solution  were  ~4  m  s-1  (Mahmud  et  al.  2011). 

2.2.  IRTF 

The  majority  of  our  NIR  observations  were  obtained  at 
the  3  m  NASA  IRTF  using  CSHELL  (Tokunaga  et  al.  1990; 
Greene  et  al.  1993).  CSHELL  is  a  long-slit  echelle  spectrograph 
(1.08-5.5 /xm)  that  uses  a  Circular  Variable  Filter  (CVF)  to 
isolate  a  single  order  onto  a  256  x  256  InSb  detector  array.  We 
used  the  CVF  to  isolate  a  50  A  segment  of  spectrum  centered 
at  2.298  /xm  (order  25).  This  region  contains  numerous  deep 
photospheric  absorption  lines  from  the  2.293  /xm  CO  v  =  2-0 
band  head  as  well  a  rich  set  of  predominately  CH4  telluric 
absorption  features  which  we  use  as  a  wavelength  reference. 
The  075  slit  yielded  a  typical  FWHM  of  2.6  pixels  (0.5  A, 
6.5  km  s“\  measured  from  arc  lamp  spectra)  corresponding 
to  a  spectral  resolving  power  of  R  ~  46,000. 

We  obtained  data  on  56  nights  between  2008  February 
and  2012  February.  These  observations  included  several  RV 
standards,  a  known  exoplanet  host,  and  nine  T  Tauri  planet 
host  candidates  (Table  1).  At  the  beginning  of  each  night,  we 
imaged  20  flat  fields,  each  with  a  20  s  integration  time,  using  a 


Table  1 

Stellar  Parameters  of  Young  Planet  Host  Candidates 
and  Spotted  Star  Standards 


Name 

V3 

(mag) 

K* 

(mag) 

Uffb 

(K) 

v  sin  *'c 
(km  s”1) 

Typeb 

(C/W)d 

Ageb 

(Myr) 

BP  Tau 

12.3 

7.7 

4060 

10.9 

C 

1.8 

Cl  Tau 

13.0 

7.8 

4060 

10.6 

C 

2.1 

DK  Tau 

12.6 

7.1 

4060 

11.5 

c 

1.0 

DN  Tau 

12.5 

8.0 

3850 

8.7 

c 

0.9 

Hubble  1 4 

12.0 

7.3 

4060 

12.9 

w 

0.7 

IQ  Tau 

14.5 

7.8 

3780 

12.0 

c 

1.4 

IW  Tau 

12.4 

8.3 

4060 

6.9 

w 

2.1 

V827  Tau 

12.2 

8.2 

4060 

19.3 

w 

2.5 

V836  Tau 

13.1 

8.6 

4060 

15.0 

c 

7.9 

Notes. 
a  SIMBAD. 

b  Palla  &  Stahler  (2002). 
c  Glebocki  &  Gnacinski  (2005). 
d  C:  classical;  W:  weak-lined. 

continuum  lamp  to  illuminate  the  entire  slit.  We  also  imaged  the 
same  number  of  20  s  dark  frames.  Additionally,  we  imaged  six 
Ar-Kr-Xe  emission  lines,  changing  the  CVF  while  maintaining 
the  grating  position,  to  determine  the  wavelength  reference.  All 
of  our  target  data  were  obtained  using  10"  nodded  pairs  to 
enable  subtraction  of  sky  emission,  dark  current,  and  detector 
bias.  Integration  times  for  each  nod  were  typically  600  s;  for 
fainter  targets,  we  took  multiple  contiguous  nod  pairs  (Crockett 
etal.  2011). 

2.3.  Keck 

We  supplemented  our  CSHELL  data  with  six  nights  of  ob¬ 
servations  using  NIRSPEC  (McLean  et  al.  1998,  2000)  between 
2007  March  and  2012  January.  NIRSPEC  is  a  NIR,  vacuum- 
cryogenic,  high-resolution,  cross-dispersed,  echelle  spectro¬ 
graph  operating  at  the  Nasymth  focus  on  the  10  m  Keck  II 
telescope.  We  used  the  N7  filter  (1.839-2.630 /xm)  with  the 
echelle  and  cross-disperser  angles  set  to  62772  and  36724,  re¬ 
spectively,  to  image  orders  30  through  35  onto  the  1024  x 
1024  InSb  detector.  This  provided  a  wavelength  coverage  of 
roughly  2.157-2.550  /xm.  The  07288  x  24"  slit  yielded  a  me¬ 
dian  FWHM  of  2.25  pixels  (0.74  A,  9.6  km  s-1,  measured  from 
arc  lamp  spectra)  corresponding  to  a  spectral  resolving  power  of 
R  ~  31,000.  To  optimize  simultaneous  telluric  and  photosphere 
line  coverage  (Section  3.3),  we  limited  subsequent  analyses  to 
order  33  (2.286-2.320  /xm). 

We  observed  RV  standards,  telluric  standards,  and  over  a 
dozen  potential  planet  hosts  in  Taurus.  At  the  beginning  of  each 
night,  we  imaged  approximately  10  flat  fields  (3  co-adds  x 
30  s)  using  a  continuum  lamp  to  illuminate  the  entire  slit.  We 
also  imaged  the  same  number  of  dark  frames  (3  co-adds  x  30  s). 
Additionally,  we  imaged  a  Ne-Ar-Xe-Kr  calibration  lamp  to 
provide  an  initial  rough  estimate  of  the  wavelength  reference. 
All  of  our  target  data  were  obtained  using  an  ABBA  10"  nod 
sequence  to  enable  subtraction  of  sky  emission,  dark  current, 
and  detector  bias.  Integration  times  were  typically  30-60  s  with 
2-3  co-adds. 

3.  DATA  REDUCTION 
3.1.  Visible  Light  Radial  Velocities 

All  spectra  were  bias-subtracted,  flat-fielded,  and  optimally 
extracted  using  an  IDL  echelle  reduction  code  (Hinkle  et  al. 


3 


The  Astrophysical  Journal,  761:164  (17pp),  2012  December  20 


Crockett  et  al. 


Table  2 

Optical  RV-Bisector  Correlations 


Target 

Pearson  Coefficient 

BP  Tau 

0.20 

Cl  Tau 

0.74 

DK  Tau 

0.60 

DN  Tau 

0.33 

Hubble  4 

0.80 

IQ  Tau 

0.28 

IW  Tau 

-0.16 

V827  Tau 

-0.88 

V836  Tau 

-0.30 

2000).  We  determined  visible  light  RVs  with  a  cross-correlation 
analysis  of  eight  echelle  orders,  each  covering  about  100  A.  The 
orders  were  chosen  for  their  high  S /N,  lack  of  stellar  emission 
lines,  and  lack  of  telluric  absorption  lines.  For  each  target’s 
cross-correlation  template  spectrum,  we  chose  the  highest  S /N 
observation.  Using  the  target  itself  as  a  template  reduces  addi¬ 
tional  uncertainty  from  any  spectral  type  mismatching.  Since 
each  order  was  cross-correlated  against  the  corresponding  order 
in  the  template,  our  RVs  are  measured  relative  to  one  obser¬ 
vation  epoch.  For  each  observation,  the  RV  and  its  associated 
uncertainty  were  taken  to  be  the  mean  and  standard  deviations, 
respectively,  of  the  RVs  from  all  utilized  echelle  orders. 

To  determine  our  long-term  RV  uncertainty,  we  observed 
several  stars  known  to  be  stable  at  a  few  m  s_1:  107  Psc, 
HD  4628,  r  Ceti,  HD  65277,  HD  80367,  and  HD  8837 1  (Nidever 
et  al.  2002;  Butler  et  al.  1996;  Cumming  et  al.  1999).  The 
standard  deviation  of  these  RV  standards  over  eight  years  is 
~140  ms-1.  We  adopt  this  value  as  our  long-term  systematic 
uncertainty.  The  final  uncertainty  for  each  observation  is  then 
140  m  s'1  added  in  quadrature  with  the  order-to-order  scatter 
(Mahmud  et  al.  2011). 

3.2.  Visible  Light  Bisector  Analysis 

We  analyzed  the  visible  light  bisectors  to  determine  the  origin 
of  each  target’s  RV  variability  (Huerta  et  al.  2008;  Mahmud  et  al. 
2011).  We  used  the  cross-correlation  function  from  each  order 
(Section  3.1)  to  determine  a  bisector  span  (the  inverse  of  the 
mean  slope  of  the  bisector);  the  mean  across  all  orders  provided 
a  mean  bisector  span  for  each  observation.  We  then  computed 
the  Pearson  linear  correlation  coefficient,  r,  between  the  bisector 
spans  and  RVs  for  each  target  (Figure  1  and  Table  2).  Those 
targets  with  no  correlation,  except  a  few  “spotted  standard  stars” 
(Section  4.2),  were  selected  for  NIR  follow-up  observations  on 
CSHELL. 

3.3.  NIR  Radial  Velocities 

All  CSHELL  and  N1RSPEC  observations  were  processed 
using  the  same  pipeline  (Crockett  et  al.  2011).  The  data  were 
dark-subtracted,  flat-fielded,  and  optimally  extracted  using  an 
IDL  implementation  of  Horne  (1986).  We  calculated  the  RVs 
using  a  spectral  modeling  technique.  This  technique  uses  two 
high-resolution  template  spectra  to  model  the  stellar  spectrum 
and  the  telluric  features.  For  the  stellar  spectrum,  we  generated 
NextGen  stellar  atmosphere  models  (Hauschildt  et  al.  1999) 
tailored  to  the  Teff,  log  g,  and  metallicity  of  our  targets.  We  then 
used  SYNTHMAG  (Piskunov  1999)  to  create  template  spectra 
from  the  NextGen  models  along  with  input  from  atomic  (Kupka 
et  al.  2000)  and  CO  (Goorvitch  1994)  line  lists. 


Table  3 

RV  Standards 


Name 

/^obs 

(JV 

(m  s-1) 

CSHELL 

GJ  281 

48 

66 

HD  65277 

9 

36 

HD  219538 

13 

55 

HD  225261 

14 

61 

NIRSPEC 

GJ  281 

9 

30 

HD  10476 

5 

129 

The  telluric  absorption  features  in  the  K  band  provide  an 
absolute  wavelength  and  instrumental  profile  reference,  similar 
in  concept  to  the  iodine  gas  cell  technique  used  in  high-precision 
optical  RV  exoplanet  surveys  (Butler  et  al.  1996).  Using  the 
atmosphere  as  a  “gas  cell”  lets  us  superimpose  a  relatively 
stable  wavelength  reference,  which  follows  the  same  optical 
path  as  the  light  from  the  science  target,  onto  our  spectra. 
This  helps  alleviate  uncertainties  introduced  by  variable  slit 
illumination,  changing  optical  path  lengths,  etc.  We  modeled  the 
telluric  features  using  the  NOAO  telluric  absorption  spectrum 
(Livingston  &  Wallace  1991). 

In  order  to  match  the  model  to  our  observations,  we  ap¬ 
plied  a  number  of  transformations  to  our  templates:  a  ve¬ 
locity  shift  and  a  power-law  scaling  factor,  rotational  and 
instrumental  broadening,  a  second-order  continuum  normal¬ 
ization,  and  a  second-order  wavelength  dispersion.  The  com¬ 
posite  model  was  then  binned  down  to  the  resolution  of  our 
CSHELL  (or  NIRSPEC)  data.  We  used  an  IDL  implementa¬ 
tion  of  the  Levenberg-Marquardt  nonlinear  least-squares  fitting 
algorithm10  to  determine  the  values  of  the  model  parameters 
which  best  reproduced  our  observations. 

We  used  Monte  Carlo  techniques  to  estimate  the  errors  in 
the  model  parameters.  For  each  observation,  we  generated  100 
simulated  observations  based  on  the  measured  noise.  We  fit  a 
model  to  each  of  the  simulated  observations  which  provided 
us  with  100  sets  of  parameters.  We  then  calculated  the  mean 
and  standard  deviation  of  the  100  results  for  each  parameter 
which  we  took  to  be  the  final  best-fit  value  and  uncertainty, 
respectively. 

The  RV  for  each  nod  position  is  the  velocity  shift  of  the  stellar 
template  relative  to  the  telluric-based  wavelength  solution.  A 
final,  nightly  RV  was  determined  by  calculating  the  average  of 
the  individual  nod  RVs,  weighted  by  their  uncertainties.  The 
final  uncertainty  in  the  nightly  RV  was  computed  by  taking  the 
weighted  standard  deviation  of  the  nod  RV s  and  dividing  by  the 
square  root  of  the  number  of  nods. 

To  assess  the  long-term  uncertainties  in  our  observations,  we 
routinely  observed  stars  known  to  have  stable  RVs  (i.e.,  er„  < 
50  m  s_  1 ).  We  observed  one  RV  standard,  GJ  28 1 ,  on  every  avail¬ 
able  night  to  ensure  that  we  had  a  consistent  test  spanning  the 
entire  project.  When  time  allowed,  we  observed  additional  RV 
standards  (HD  65277,  HD  219538,  HD  225261).  Table  3  sum¬ 
marizes  the  rms  RV  s  for  each  standard.  We  use  the  rms  of  GJ  28 1 , 
66  m  s-1,  as  the  canonical  value  of  our  long-term  CSHELL 
stability  since  it  is  the  only  target  observed  consistently  on 
all  runs. 


1°  mpfit,  written  by  Craig  Markwardt 
(http://cow.physics.wisc.edu/~craigm/idl/fitting.html) 
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Figure  1.  Optical  RV  vs.  bisector  span  for  Taurus  targets  with  correlation  coefficients.  All  data  were  acquired  on  the  McDonald  2.7  m  telescope.  Strong  correlations 
suggest  that  the  optical  RV  modulation  is  spot  driven. 


GJ  281  was  also  observed  on  most  nights  with  NIRSPEC, 
in  addition  to  HD  10476  (Table  3).  We  use  the  rms  of  GJ  281. 
30  m  s_I,  to  quantify  our  night-to-night  uncertainties  because, 
as  with  CSHELL,  it  is  the  most  frequently  observed  standard 
with  the  longest  time  baseline.  Using  GJ  281  to  assess  our 
NIRSPEC  uncertainties  also  provides  a  common  benchmark 
for  comparisons  to  our  CSHELL  observations. 

4.  RV  ANALYSIS 
4.1.  Gl  86 

In  Crockett  et  al.  (2011),  we  demonstrated  the  ability  of 
CSHELL  to  detect  giant  planets  by  confirming  the  exoplanet 
companion  around  Gl  86  (Queloz  et  al.  2000;  Figueira  et  al. 
2010).  Figure  2  plots  new  CSHELL  RVs  with  those  from 
Queloz  et  al.  and  Figueira  et  al.  along  with  an  orbital  fit  to 
all  data  sets.  Our  error  bars  are  determined  by  adding  the 
uncertainties  from  our  RV  algorithm  in  quadrature  with  the 
66  m  s-1  systematic  uncertainty  (Section  3.3).  Our  best-fit  orbit 
(Table  4)  is  very  close  to  the  Queloz  et  al.  solution.  Our  data  show 
good  agreement  with  previously  published  values  exhibiting  a 
~36  m  s'1  standard  deviation  in  the  O—C  residuals. 

The  CSHELL  RVs  in  Figure  2  have  been  corrected  for  a 
long-term  linear  drift  in  this  object,  first  reported  in  Queloz 


Figure  2.  RVs  for  known  exoplanet  Gl  86  b.  Filled  red  circles  are  CORALIE 
RVs  from  Queloz  et  al.  (2000),  blue  stars  are  CRIRES  RVs  from  Figueira 
et  al.  (2010),  and  open  green  squares  are  our  K- band  CSF1ELL  RVs  from  2008 
November  to  2011  November.  The  CSHELL  RVs  have  been  corrected  for  a 
long-term  linear  drift  (see  Figure  3).  The  black  line  is  an  orbit  model  fit  to  the 
combined  data  sets.  The  standard  deviation  of  our  residuals  is  36  m  s_1 . 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Table  4 

G1  86  b  Parameters 


Parameter 

Value 

Units 

K 

381.4 

m  s_1 

P 

15.764 

days 

e 

0.05 

To 

2451146.836 

JD 

CO 

273.4 

deg 

dv 

-0.22 

m  s”1  d”1 

o  - 

"  ^Queloz 

7.9 

m  s_1 

o  - 

-  Cpigueria 

24.2 

m  s_1 

o  - 

"  ^Crockett 

36.1 

m  s_1 

et  al.  (2000).  Figure  3  plots 
subtracted  and  a  linear  fit 

our  RV s  with  the  planet  modulation 
to  the  residuals.  A  Monte  Carlo 

analysis  shows  that  our  data  are  consistent  with  a  0.22  ± 
0.04  m  s”1  day-1  drift.  This  is  shallower  than  the  0.36  m  s'1 
day”1  drift  reported  by  Queloz  et  al.  Both  data  sets  suggest 
another  companion  with  an  orbital  period  on  the  order  of 
hundreds  of  years  (a  >  20  AU).  Our  results  demonstrate  the 
ability  of  CSHELL  to  detect  not  only  short-period  Jupiter-mass 
planets  but  also  secular  trends  from  more  distant  companions. 

4.2.  Discussion  of  Individual  Taurus  Targets 
4.2.1.  Oven’iew 

From  the  56  stars  on  our  target  list  with  a  complete  set  of 
12  visible  light  observations,  only  13  have  met  our  criteria 
for  NIR  follow-up  observations.  Eight  of  those  currently  have 
enough  high  S /N  observations  to  make  meaningful  comparisons 
between  the  two  wavelength  bands  (the  remaining  five  will 
continue  to  be  monitored  in  upcoming  observing  runs).  In 
addition  to  these  “planet  host  candidates,”  we  observed  a  known 
spotted  star,  V827  Tau,  to  help  us  understand  how  spots  will 
affect  the  long  wavelength  RV  behavior.  These  targets,  along 
with  relevant  stellar  parameters,  are  listed  in  Table  1.  The 
remainder  of  the  sample  will  be  discussed  in  another  paper. 

For  each  of  our  targets,  we  calculated  the  peak-to-peak  am¬ 
plitude  of  the  RV  variation  in  both  optical  and  NIR  wavelengths, 
looked  for  periodicity  at  each  wavelength,  and  fitted  Keplerian 
models  to  significant  periods  in  both  bands.  To  estimate  the 
amplitudes  and  their  associated  uncertainties,  we  generated  105 
data  sets  for  every  target,  replacing  each  RV  with  a  Gaussian 
random  number;  the  mean  and  standard  deviation  were  set  to  the 
observed  RV  and  its  associated  error  bar,  respectively.  For  each 
simulated  data  set,  we  recorded  the  peak-to-peak  amplitudes 
for  the  optical  and  NIR  RVs  ( A0  and  A/,  respectively)  and  the 
amplitude  ratio  between  the  two  wavelengths  <R,\)'.  we  did  the 
same  for  the  RV  standard  deviations  in  the  two  spectral  bands. 
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Figure  3.  G1  86  CSHELL  RVs  after  subtracting  the  planet-induced  modulation. 
The  dashed  line  is  a  linear  fit  to  the  residuals  of  —0.22  m  s”1  day”1 .  A  similar 
long-term  trend  was  presented  by  Queloz  et  al.  (2000)  and  indicates  the  presence 
of  a  second,  more  distant,  companion. 


We  then  fit  Gaussian  functions  to  histograms  of  all  measured 
values.  The  mean  of  each  Gaussian  was  taken  to  be  the  canoni¬ 
cal  value  of  that  parameter  with  error  bars  equal  to  the  standard 
deviation  (Table  5  and  Figure  4). 

To  look  for  periodicity  in  our  RV  measurements,  we  used 
an  IDL  implementation  of  the  CLEAN  algorithm  (Roberts 
et  al.  1987).  CLEAN  attempts  to  remove  the  couplings  between 
physical  periods  and  their  aliases  by  deconvolving  the  spectral 
window  function  from  the  discrete  Fourier  transform.  We 
calculated  the  power  spectrum  for  each  target  using  the  RVs 
from  the  (1)  optical  observations  only,  (2)  NIR  observations 
only,  and  (3)  combined  optical  and  NIR  data  sets. 

We  also  fit  a  Keplerian  model  to  both  the  optical  and  NIR 
data  for  each  target  phased  to  significant  periods  found  in 
the  periodicity  analysis.  To  help  us  quantify  differences  in  the 
optical  and  NIR  RVs,  we  started  by  fitting  five  free  parameters 
( K ,  e,  Tq,  a>,  and  Vo)  to  the  optical  data.  The  period  was  fixed  to 
the  strongest  peak  in  the  optical  power  spectrum.  Since  we  were 
primarily  interested  in  how  the  RV  amplitude  changes  at  the 
longer  wavelengths,  we  repeated  the  fit  to  the  NIR  data  while 
allowing  only  the  amplitude  to  vary.  All  of  the  other  parameters 
were  fixed  at  the  best-fit  values  for  the  optical  RVs.  For  targets 
that  showed  strong  periodicity  in  the  NIR,  we  performed  the 
additional  step  of  reversing  this  procedure:  we  fit  a  model  to  the 
NIR  RVs  with  the  period  fixed  to  the  peak  of  the  NIR  power 
spectrum.  We  then  fit  that  model  to  the  optical  data,  once  again 
allowing  only  the  amplitude  to  vary.  For  clarity,  we  refer  to 


Table  5 

Optical  and  NIR  RV  Amplitudes,  Standard  Deviation,  and  Ratios  Between  Wavelengths 


Target 

Amplitudes 

Standard  Deviations 

Optical 

NIR 

Ratio 

Optical 

NIR 

Ratio 

BP  Tau 

2.68 

± 

0.26 

1.70 

± 

0.26 

1.45 

± 

0.28 

0.64 

± 

0.03 

0.43 

± 

0.04 

1.43 

± 

0.15 

Cl  Tau 

2.80 

± 

0.30 

1.90 

± 

0.17 

1.41 

± 

0.20 

0.95 

± 

0.08 

0.63 

± 

0.05 

1.50 

± 

0.17 

DK  Tau 

5.64 

± 

0.31 

3.46 

± 

0.24 

1.50 

± 

0.16 

1.58 

± 

0.07 

0.85 

± 

0.05 

1.80 

± 

0.13 

DN  Tau 

2.12 

± 

0.26 

0.87 

± 

0.13 

2.27 

± 

0.42 

0.59 

± 

0.05 

0.27 

± 

0.03 

2.13 

± 

0.30 

Hubble  I  4 

3.34 

± 

0.22 

1.51 

± 

0.14 

2.11 

± 

0.23 

1.07 

± 

0.04 

0.41 

± 

0.03 

2.54 

± 

0.19 

IQ  Tau 

4.97 

it 

0.35 

2.59 

± 

0.54 

1.75 

± 

0.41 

1.19 

± 

0.07 

0.77 

± 

0.11 

1.48 

± 

0.23 

IW  Tau 

1.40 

± 

0.18 

0.87 

± 

0.20 

1.44 

± 

0.37 

0.36 

± 

0.03 

0.25 

± 

0.04 

1.36 

± 

0.23 

V827  Tau 

5.40 

± 

0.86 

1.16 

± 

0.31 

4.24 

± 

1.13 

1.82 

± 

0.06 

0.37 

± 

0.09 

4.63 

± 

1.02 

V836  Tau 

2.90 

± 

0.42 

1.06 

± 

0.30 

2.45 

± 

0.74 

0.80 

± 

0.08 

0.35 

± 

0.09 

2.11 

± 

0.54 
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Figure  4.  Optical  and  NIR  RV  amplitudes  (top)  and  standard  deviations  (bottom).  The  dashed  lines  in  each  plot  show  optical-NIR  ratios  of  one,  two,  and  four.  Targets 
with  primarily  planet-induced  modulation  should  fall  along  the  line  A  =  1,  whereas  those  with  larger  A  values  have  stronger  spot  signatures. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


these  four  models  as  follows:  optical  Kepler  fit  to  the  optical 
data  (OKOD),  optical  Kepler  fit  to  the  infrared  data  (OKID), 
infrared  Kepler  fit  to  the  infrared  data  (IKID),  and  infrared 
Kepler  fit  to  the  optical  data  (IKOD). 

4.2.2.  BP  Tau 

BP  Tau  is  our  most  intensively  observed  target.  Initial  NIR 
observations  strongly  supported  a  giant  planet  hypothesis; 
continued  observing  has  failed  to  uphold  that.  BP  Tau  exhibits 
significant  optical  RV  variability  with  Ao  =  2.7  km  s-1  and  a 
period  of  8.28  days;  the  optical-NIR  peak-to-peak  RV  amplitude 
ratio  (Ra)  is  1.45.  The  Kepler  fits  to  the  phase-folded  RV 
curves  (OKOD  and  OKID,  Figure  5)  suggest  a  spot-driven 
interpretation  of  the  RV  modulation.  The  ratio  of  optical-NIR 
Kepler  model  amplitudes  (Rk)  is  ~4.5.  The  NIRSPEC  data,  in 
particular,  show  only  slight  variation  at  about  the  2cr  level. 

The  optical  RV  period  closely  matches  photometric  periods 
reported  in  the  literature  of  8.3  and  8.19  days  (Richter  et  al. 
1992;  Percy  et  al.  2006).  High  time-resolution  photometric  data 
from  the  Trans- Atlantic  Exoplanet  Survey  (TrES)  also  show  a 
comparable  period  of  8.309  days  (Xiao  et  al.  2012).  However, 
periods  of  6.1, 7.6,  and  7.7  days  have  also  been  identified  (Simon 
et  al.  1990;  Vrba  et  al.  1986;  Osterloh  et  al.  1996,  respectively). 


The  correlation  with  known  photometric  periods,  the  ratio  of 
both  peak-to-peak  and  model  amplitudes,  as  well  as  the  lack  of 
any  other  significant  periodicity  in  the  NIR  RVs,  all  point  to  BP 
Tau’s  observed,  short-period  RV  variation  being  spot  induced. 

4.2.3.  Cl  Tau 

Cl  Tau  exhibits  strong  optical  RV  periodicity  at  6.47  days. 
The  NIR  data  do  not  phase  meaningfully  to  this  period  but  do 
exhibit  significant  variability  (Ra  —  1.41;  Figure  6).  In  the 
combined  CSHELL  and  NIRSPEC  data,  we  find  a  peak  in  the 
periodogram  at  10.87  days  which  is  not  seen  in  the  optical  data. 
Figure  7  presents  the  same  optical  and  NIR  RVs  phased  to  this 
other  period.  In  this  case,  we  freely  fit  all  parameters  in  the  NIR 
Kepler  fit  and  then  fit  the  optical  data  by  fixing  all  parameters 
except  the  amplitude.  The  optical  fit  is  clearly  nonsensical;  there 
is  no  periodicity  at  10.87  days  in  the  optical  RVs. 

Assuming  the  IKID  models  a  true  companion,  we  subtracted 
the  planet  modulation  from  both  the  optical  and  NIR  RVs. 
Figure  8  illustrates  the  modified  RVs  phased  to  the  6.47  day 
period.  The  optical  RV s  demonstrate  significant  improvement.  A 
Keplerian  model  to  the  new  optical  RVs  also  shows  an  excellent 
fit.  In  this  scenario,  the  optical  Keplerian  parameters  are  most 
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BP  Tau,  P  =  8.28  days 


Figure  5.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  BP  Tau  phased  to  8.28  days.  On  the  right,  red  circles  are  CSHELL  data  and  green  squares  are  NIRSPEC.  The 
black  line  is  a  best-fit  Kepler  model.  The  large  change  in  amplitude  between  the  optical  and  NIR  Kepler  models  indicates  that  the  RV  modulation  is  most  likely  spot 
induced. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Cl  Tau,  P  =  6.47  days 


Figure  6.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  Cl  Tau  phased  to  6.47  days.  On  the  right,  red  circles  are  CSHELL  data  and  green  squares  are  NIRSPEC.  The 
black  line  is  a  best-fit  Kepler  model.  The  extremely  poor  fit  to  the  NIR  data  indicates  that  the  periodicity  seen  in  the  optical  RVs  is  not  evident  in  the  NIR.  The  NIR 
and  optical  RVs  are  most  likely  caused  by  different  physical  mechanisms  (see  Figure  7). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


likely  modeling  the  spot  modulation.  The  same  model  fit  to  the 
subtracted  NIR  RVs,  allowing  only  the  amplitude  to  vary,  is  flat. 

Although  for  Cl  Tau  we  have  fewer  than  12  optical  observa¬ 
tions,  we  see  a  strong  correlation  ( r  —  0.74)  between  the  optical 
RVs  and  bisector  spans.  Furthermore,  we  observed  significant 
variability  in  the  K  band  with  a  different  periodicity  than  the 
optical  RVs.  These  data  suggest  that  a  heavily  spotted  star  with 
a  planet  may  exhibit  optical  variability  that  is  entirely  consis¬ 
tent  with  a  spot-only  model  for  the  RV  modulation.  Large  spots 
may  mask  planet-induced  Doppler  shifts.  Modeling  is  needed 
to  verify  this  possibility  (C.  J.  Crockett  et  al.  in  preparation). 


4.2.4.  DK  Tau 

DK  Tau  exhibits  large  RV  amplitudes  in  both  wavelengths: 
5.6  km  s-1  in  the  optical  and  3.5  km  s-1  in  the  NIR.  An 
amplitude  ratio  of  Ra  =  1.5  and  a  linear  correlation  coefficient 
of  r  —  0.6  are  consistent  with  spot-induced  RV  variability.  It 
is  not  clear,  however,  that  a  simple  spot  model  is  sufficient  to 
explain  the  observed  behavior.  Figure  9  shows  the  optical  and 
NIR  RVs  phase  folded  to  the  strongest  period  in  the  optical 
power  spectrum,  P  =  4.83  days.  There  is  significant  scatter 
in  the  optical  RVs.  After  subtracting  the  OKOD,  the  standard 
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Cl  Tau,  P  =  10.87  days 


Figure  7.  Same  as  Figure  6,  but  phased  to  the  strongest  NIR  period  at  10.87  days.  A  Kepler  model  is  unable  to  fit  the  optical  data  at  this  period  (black  line,  left)  but 
does  provide  a  good  fit  to  the  NIR  RVs.  This,  coupled  with  the  6.47  day  fit  and  optical  bisector  analysis,  suggests  that  the  optical  RVs  are  predominately  spot  driven, 
whereas  the  NIR  RVs  arise  from  a  different  mechanism,  possibly  a  companion. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


> 

DC 


as 

<D 

DC 


Cl  Tau,  P  =  6.47  days 


Figure  8.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  Cl  Tau  after  the  best-fit  NIR  Kepler  model  (Figure  7)  have  been  subtracted  from  both  data  sets.  On  the  right,  red 
circles  are  CSHELL  data  and  green  squares  are  NIRSPEC.  The  black  line  is  a  new  Kepler  model  fit  to  the  modified  data.  The  optical  RV  modulation  is  much  cleaner 
than  Figure  6  and  the  NIR  RV  residuals,  as  expected,  show  little  variation.  This  supports  our  suggestion  that  the  optical  RV  modulation  is  the  superposition  of  both 
spot-  and  companion-induced  variabilities. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


deviation  in  the  residuals  is  ~1  km  s~  1 .  Furthermore,  the  NIR 
RV s  do  not  phase  to  this  period  at  all;  the  amplitude  of  the  OKID 
is  flat  with  ~1  km  s-1  residuals. 

The  strongest  period  in  the  NIR  RVs  is,  as  for  Cl  Tau 
(Section  4.2.3),  significantly  different  than  that  seen  in  the 
optical:  14.93  days.  Figure  10  shows  RV  data  from  both 
wavelengths  phased  to  this  period.  The  IKID  is  a  better  fit 
(cr o-c  ~  760  m  s”1)  than  the  OKID,  but  the  NIR  RVs  still 
show  considerable  scatter.  The  optical  data  do  not  phase  to  this 
period. 


4.2.5.  DN  Tau 

DN  Tau  was  interpreted  as  a  spotted  star  in  Prato  et  al.  (2008) 
on  the  basis  of  eight  NIR  RVs  obtained  in  2008  February. 
Revisiting  these  data  with  our  improved  NIR  reduction  pipeline 
hinted  at  a  possibly  more  complex  system  so  DN  Tau  was  added 
back  to  the  CSHELL  target  list  for  more  intensive  follow¬ 
up.  DN  Tau  shows  little  RV-bisector  correlation  ( r  —  0.33) 
and  considerable  scatter  in  the  NIR  (A/  ss  870  m  s-1).  The 
amplitude  ratio  ( Ra  ~  2.3),  however,  seems  to  again  place  this 
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DK  Tau,  P  =  4.83  days 


Figure  9.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  DK  Tau  phased  to  4.83  days.  On  the  right,  red  circles  are  CSHELL  data  and  green  squares  are  NIRSPEC.  The 
black  line  is  a  best-fit  Kepler  model.  Despite  being  the  strongest  period  in  the  optical  periodogram,  the  fit  is  poor  at  both  wavelengths. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


DK  Tau,  P  =  14.93  days 


Figure  10.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  DK  Tau  phased  to  14.93  days.  On  the  right,  red  circles  are  CSHELL  data  and  green  squares  are  NIRSPEC.  The 
black  line  is  a  best-fit  Kepler  model.  The  optical  RVs  show  no  periodicity  with  this  phasing  (as  evidenced  by  the  attempt  at  Kepler  fit).  The  NIR  RV  fit  is  much  better 
than  the  one  seen  at  4.83  days,  though  with  considerable  scatter.  This  may  be  indicative,  like  Cl  Tau,  of  different  physical  mechanisms  driving  the  optical  and  NIR  RV 
modulation. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


target  firmly  in  the  spot-induced  category.  Figure  1 1  shows  the 
optical  and  NIR  data  both  phased  to  6.4  days  (the  strongest  peak 
in  the  optical  periodogram)  along  with  the  OKOD  and  OKID. 
The  ratio  of  Keplerian  model  amplitudes  ( R K  2.7)  further 
supports  the  spot  hypothesis. 

However,  interpretation  of  DN  Tau,  similar  to  Cl  Tau  (Sec¬ 
tion  4.2.3),  is  muddied  by  a  second  significant  period  in  the 
NIR  data  at  2.65  days  (Figure  12).  To  fit  the  NIR  data,  we  fol¬ 
lowed  the  same  approach  taken  for  Cl  Tau.  The  IKOD  confirms 
the  optical  power  spectrum  analysis:  there  is  no  periodicity  at 


2.65  days  in  the  optical  RVs.  This  is  another  target  that  will 
benefit  from  increased  NIR  scrutiny. 

4.2.6.  Hubble  1 4 

Mahmud  et  al.  (2011)  reported  that  the  RV  modulation  of 
Hubble  I  4  is  consistent  with  a  large  polar  spot  on  a  nearly 
pole-on  star.  Here,  we  refine  that  analysis  with  additional 
NIRSPEC  RVs  and  improvements  to  the  NIR  spectrum  model¬ 
ing  (Figure  13). 
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DN  Tau,  P  =  6.40  days 


Figure  11.  Optical  RVs  (left)  and  CSHELL  RVs  (right)  of  DN  Tau  phased  to  6.40  days.  The  black  line  is  a  best-fit  Kepler  model.  The  reduced  amplitude  of  the  NIR 
fit  indicates  that  this  RV  variability  is  primarily  spot  induced. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


DN  Tau,  P  =  2.65  days 


Figure  12.  Optical  RVs  (left)  and  CSHELL  RVs  (right)  of  DN  Tau  phased  to  2.65  days.  The  black  line  is  a  best-fit  Kepler  model.  This  period  is  the  strongest  seen  in 
the  NIR  periodogram.  The  optical  data  do  not  phase  to  this  period.  The  NIR  fit  hints  at  a  different  mechanism  driving  the  NIR  variability,  but  further  investigation  is 
needed. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


We  previously  published  a  rotation  period  of  1.5459  days, 
arguing  that  the  slightly  stronger  2.81  day  period — seen  in  a 
Scargle  periodogram  analysis — was  most  likely  an  alias  of  the 
true  period.  We  based  this  argument  on  an  observed  reduc¬ 
tion  in  the  2.81  day  power  in  the  CLEAN  periodogram  and 
on  a  published  SuperWASP  photometric  period  of  1.5483  days 
(Norton  et  al.  2007).  Xiao  et  al.  (2012)  also  report  a  photo¬ 
metric  period  of  1.547  days  seen  in  TrES  data  (private  com¬ 
munication).  The  addition  of  eight  NIRSPEC  observations  has 
strengthened  the  case  for  a  1.5459  day  period.  Revisiting  the 
periodogram  analysis  of  Mahmud  et  al.  (2011)  with  these 


new  data,  the  2.81  day  period  no  longer  appears  in  the  power 
spectrum. 

Not  only  does  Hubble  I  4  support  our  argument  that  optical 
observations  are  not  sufficient  when  looking  for  planets  around 
young,  active  stars,  it  argues  that  relying  solely  on  NIR  observa¬ 
tions  is  also  not  conclusive.  Hubble  1 4  was  the  first  target  in  our 
survey  to  exhibit  resolved  RV  modulation  in  the  K  band.  These 
data  alone  are  consistent  with  the  presence  of  a  giant  planet. 
However,  comparison  to  the  much  larger  variability  in  the  opti¬ 
cal  data,  with  the  same  periodicity  and  phase,  demonstrates  that 
the  RV  modulation  of  Hubble  I  4  is  activity  induced. 
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Hubble  4,  P  =  1 .55  days 


Figure  13.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  Hubble  I  4  phased  to  1.55  days.  On  the  right,  red  circles  are  CSHELL  data  and  green  squares  are  NIRSPEC. 
The  black  line  is  a  best-fit  Kepler  model.  The  sharply  reduced  amplitude  in  the  NIR  and  consistent  phasing  between  the  two  wavelengths  is  a  strong  argument  for 
spot-driven  RV  modulation. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


IQ  Tau,  P  =  1.27  days 


Figure  14.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  IQ  Tau  phased  to  1.27  days.  On  the  right,  red  circles  are  CSHELL  data  and  green  squares  are  NIRSPEC.  The 
black  line  is  a  best-fit  Kepler  model.  The  NIR  RVs  show  little  variability  outside  their  error  bars,  indicating  that  the  optical  RV  variability  is  most  likely  spot  induced. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


4.2.7.  IQ  Tau 

At  a  typical  V  magnitude  of  14.5,  IQ  Tau  is  one  of  our 
faintest  McDonald  targets.  The  optical  RVs  show  periodicity 
at  1.27  days  with  Ap  5  km  s-1.  The  NIR  RVs  do  not 
phase  to  this  period  and  have  a  lower  peak-to-peak  amplitude  of 
2.6  km  s-1  (Figure  14).  The  reduction  in  RV  variability  at  longer 
wavelengths  is  strong  evidence  that  the  optical  RV  behavior  is 
spot  induced.  The  NIR  power  spectrum  does  have  a  peak  at 
9.5  days,  shown  in  Figure  15.  However,  the  IKID  seems  to  be 
biased  toward  a  slight  dip  from  one  data  point  at  phase  0. 1 .  The 
NIRSPEC  RVs  are  not  consistent  with  this  drop  and,  in  fact, 
show  very  little  variation  overall. 


During  our  observing  runs  at  McDonald  in  early  2012,  imag¬ 
ing  at  the  0.8  m  revealed  dramatic  variations  in  magnitude.  Four 
months  of  intensive  optical  photometry  from  McDonald  Ob¬ 
servatory,  Lowell  Observatory,  and  the  U.S.  Naval  Observatory 
Flagstaff  Station  (to  be  published  in  a  later  paper)  show  extreme 
periodic  (~6  days)  BVRI  magnitude  variations  up  to  2.5  mag. 
We  have  not  yet  determined  any  reasonable  explanation  for  the 
photometric  variability  of  this  system. 

4.2.8.  IW  Tau 

IW  Tau  is  another  target  that  has  comparable  RV  amplitudes 
in  optical  and  NIR  wavelengths  ( R A  ~  1 .4)  and  shows  evidence 
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IQ  Tau,  P  =  9.50  days 


Figure  15.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  IQ  Tau  phased  to  9.50  days.  On  the  right,  red  circles  are  CSHELL  data  and  green  squares  are  NIRSPEC.  The 
black  line  is  a  best-fit  Kepler  model.  The  9.5  day  period  is  the  strongest  seen  in  the  NIR  periodogram.  The  optical  RVs  show  considerable  scatter  around  the  model, 
most  likely  caused  by  star  spots.  The  NIR  fit  is  driven  mostly  by  the  CSHELL  RVs;  the  NIRSPEC  RVs  do  not  show  the  same  degree  of  variability. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


IW  Tau,  P  =  2.95  days 


Relative  Phase 

Figure  16.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  IW  Tau  phased  to  2.95  days.  On  the  right,  red  circles  are  CSHELL  data,  green  squares  are  NIRSPEC.  The  black 
line  is  a  best-fit  Kepler  model.  The  extremely  poor  fit  on  the  right  shows  that  the  strongest  optical  period  does  not  show  itself  in  the  NIR  data.  There  may  be  a  different 
physical  mechanism  behind  the  large  variability  seen  in  the  two  wavelengths. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

of  different  periods  in  the  two  bands.  At  the  strongest  optical 
period  (2.95  days),  the  optical  RVs  fit  reasonably  well  to  an 
eccentric  Kepler  model,  whereas  the  NIR  RVs  do  not  phase 
to  this  period  at  all  (Figure  16).  The  strongest  NIR  period  is 
19.57  days  (Figure  17).  Phased  to  this  period,  the  optical  RVs 
show  little  periodicity  while  the  NIR  RVs,  though  plagued  by 
significant  scatter,  show  an  adequate  fit  to  a  Keplerian  model. 

4.2.9.  V827  Tau 

We  added  V827  Tau  to  our  target  list  as  a  “spotted  standard 
star”  (Prato  et  al.  2008).  We  have  used  it  as  a  benchmark  for 
how  a  known  spotted  star  should  behave  at  both  wavelengths. 

We  see  strong  periodicity  in  the  optical  RVs  at  3.62  days.  Optical 
RVs  and  bisector  span  show  a  strong  correlation  (Figure  1)  with 
a  linear  correlation  coefficient  of  r  =  —0.88,  consistent  with 


spot-driven  RV  modulation.  Figures  4  and  18  show  a  clean 
separation  between  the  optical  and  NIR  ranges  with  Ra  ~  4.2. 
The  Keplerian  fits  in  Figure  1 8  also  show  a  stark  contrast.  The 
OKOD  has  an  amplitude  of  2.28  km  s_  1 ,  whereas  the  amplitude 
for  the  OKID  is  only  0.15  km  s-1 — a  ratio  of  ~15! 

V827  Tau  is  another  excellent  demonstration  of  the  need  for 
follow-up  of  optical  observations  at  longer  wavelengths:  the  op¬ 
tical  data  alone  fit  reasonably  well  to  a  Keplerian  model,  whereas 
the  NIR  data  show  little  variation  outside  of  the  uncertainties. 
Without  the  NIR  data  (and  bisector  analysis),  it  would  be  easy 
to  conclude  that  the  RV  variability  was  planet  induced. 

4.2.10.  V836  Tau 

V836  Tau  was  also  originally  presented  in  Prato  et  al. 
(2008)  with  the  conclusion  that,  despite  a  lack  of  RV-bisector 
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IW  Tau,  P  =  19.57  days 


Figure  17.  Optical  RVs  (left)  and  NIR  RVs  (right)  of  IW  Tau  phased  to  19.57  days.  On  the  right,  red  circles  are  CSHELL  data  and  green  squares  are  NIRSPEC.  The 
black  line  is  a  best-fit  Kepler  model.  At  this  period,  the  optical  RVs  do  not  phase  nor  do  they  fit  a  Keplerian  model.  The  NIR  RVs  exhibit  much  better  phase  coherence 
than  that  seen  at  2.95  days.  There  is,  however,  large  scatter  seen  in  the  NIR  residuals. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  18.  Optical  RVs  (left)  and  CSHELL  RVs  (right)  of  V827  Tau — our  “spotted  standard” — phased  to  3.62  days.  The  black  line  is  a  best-fit  Kepler  model.  The 
drastic  change  in  amplitude  between  the  two  wavelengths  is  the  classic  signature  of  a  heavily  spotted  star. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


V 827  Tau,  P  =  3.62  days 


correlation  ( r  —  —0.3),  the  RV  modulation  was  spot  induced  and 
not  the  result  of  an  unseen  companion.  The  numbers  presented 
here  are  updated  based  on  additional  NIR  observations  and 
improvements  to  the  CSHELL  data  reduction  pipeline.  Despite 
significant  optical  RV  variability  (A<?  =  2.9  km  s_1),  the  NIR 
scatter  is  almost  entirely  within  the  CSHELL  error  bars.  With 
Ra  ~  2.45  and  RK  «  2.7  (P  =  2.48  days;  Figure  19),  the  data 
support  the  spot  hypothesis. 

We  do  see  some  periodicity  in  the  NIR  data  at  3.55  days 
(Figure  20).  We  repeated  the  approach  of  Cl  Tau  and  DN  Tau  by 
fixing  the  optical  model  parameters  to  those  determined  by  the 
NIR  fit  (Section  4.2.5).  The  IKID  fits  the  NIR  data  better  than 


the  OKID;  the  infrared  model  does  not  fit  the  optical  RVs  at 
all.  However,  the  low  level  of  NIR  variability  prevents  us  from 
saying  anything  definitive.  This  is  yet  another  target  that  can 
greatly  benefit  from  intensive  NIR  follow-up,  particularly  with 
NIRSPEC. 

5.  DISCUSSION 

Our  initial  hypothesis  that  spot-induced  RV  modulation 
should  decrease  at  NIR  wavelengths  was  primarily  based  on 
multiwavelength  photometric  surveys  (e.g.,  Carpenter  et  al. 
2001;  Vrba  et  al.  1986).  However,  there  was  a  lack  of  long-term 
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V836  Tau,  P  =  2.48  days 


Figure  19.  Optical  RVs  (left)  and  CSHELL  RVs  (right)  of  V836  Tau  phased  to  2.48  days.  The  black  line  is  a  best-fit  Kepler  model.  There  is  little  variability  outside 
the  error  bars  in  the  NIR  RVs,  suggesting  that  this  target’s  RV  modulation  is  spot  induced. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


V836  Tau,  P  =  3.55  days 


Figure  20.  Optical  RVs  (left)  and  CSHELL  RVs  (right)  of  V836  Tau  phased  to  3.55  days.  The  black  line  is  a  best-fit  Kepler  model.  Phased  to  a  longer  period,  the 
optical  RVs  lose  phase  coherence,  whereas  the  NIR  RVs  show  a  slightly  better  fit  to  a  Kepler  model.  The  amplitude  is,  however,  at  the  same  scale  as  the  error  bars. 
The  reasonable  NIR  fit  suggests  that  higher  precision  RV  follow-up  (i.e.,  with  NIRSPEC)  is  warranted  before  the  companion  hypothesis  can  be  ruled  out. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


NIR  spectroscopic  monitoring  of  T  Tauri  stars  on  which  to 
rely.  Our  own  spot  simulations  (to  be  published  in  a  later 
paper,  but  see  also  Huerta  et  al.  2008;  Mahmud  et  al.  2011),  in 
addition  to  other  published  models  (Desort  et  al.  2007;  Reiners 
et  al.  2010;  Ma  &  Ge  2012),  predict  that  NIR  RV  monitoring 
reduces  the  stellar  jitter  of  young  stars  by  a  factor  of  two 
to  three.  Data  from  several  young  stars — TW  Hya  (Huelamo 
et  al.  2008),  V827  Tau  (Prato  et  al.  2008),  and  Hubble  I  4 
(Mahmud  et  al.  2011) — confirm  this  prediction.  Furthermore, 
observations  of  Hubble  I  4  demonstrate  another  important 
point:  NIR  observations  alone  are  not  sufficient.  Even  in  the 
K  band,  some  young  stars  exhibit  RV  amplitudes  on  the  order 


of  hundreds  of  meters  per  second  with  good  phase  coherence 
over  several  years.  Only  by  comparing  amplitudes  in  at  least 
two  widely  separated  bands  can  we  hope  to  distinguish  between 
stellar  activity  and  true  companions. 

Observations,  however,  have  forced  us  to  consider  that  a  cool 
spot  model,  at  times,  may  be  too  simplistic.  BP  Tau,  for  example, 
exhibits  high  levels  of  NIR  RV  variability  (~1.7  km  s_I  peak- 
to-peak)  which  does  not  phase  to  the  optical  data,  nor  shows 
any  strong  periodicity  of  its  own,  and  cannot  be  explained 
solely  by  the  presence  of  star  spots.  Eiroa  et  al.  (2002)  present 
results  from  simultaneous  optical  and  NIR  photometry  of  PMS 
stars  which  suggest  alternate  explanations.  They  divide  their 
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sample  into  two  groups:  (1)  stars  with  similar  variability  trends 
in  optical  and  NIR  bands  and  (2)  stars  with  different  trends. 
Those  stars  in  (1)  can  be  explained  by  the  presence  of  star 
spots  or  variable  obscuration  (such  as  from  a  dusty  disk).  Stars 
in  (2)  are  more  challenging  and  not  all  show  the  expected 
decrease  in  amplitude  variation  with  wavelength.  In  some 
cases,  the  NIR  variability  is  greater  than  that  in  the  optical; 
other  targets  get  brighter  in  one  band  while  getting  dimmer  in 
another.  This  suggests  that  for  some  PMS  stars,  the  optical  and 
NIR  variability  have  different  underlying  physical  mechanisms. 
While  the  optical  light  is  dominated  by  the  photosphere,  the  NIR 
flux  has  some  contribution  from  the  circumstellar  disk.  They 
estimate  that  the  disk  contributes  more  than  50%  of  the  K-band 
flux  in  their  sample  and  thus  conclude  that  the  NIR  variability 
of  stars  with  a  NIR  excess  predominately  tracks  changes  in 
the  disk. 

In  the  case  of  BP  Tau,  Folha  &  Emerson  (1999)  calculate  a 
K-band  veiling  of  0.8  ±  0.3  based  on  high  resolution  spec¬ 
troscopy  near  the  2. 1661  /xm  Bry  line.  This  implies  that  roughly 
44%  of  the  A'-band  light  is  non-photospheric.  We  therefore  can¬ 
not  rule  out  the  possibility  that  our  /("-band  observations  are 
impacted  by  either  changes  in  the  disk  structure  or  shocks  from 
the  base  of  magnetospheric  accretion  columns.  Muzerolle  et  al. 
(2003)  present  models  that  support  this  hypothesis.  They  show 
that  the  NIR  excess  in  a  sample  of  classical  T  Tauri  stars  (in¬ 
cluding  BP  Tau)  is  consistent  with  a  single-temperature  black- 
body  with  T  ~  1400  K.  The  amount  of  excess  also  roughly 
correlates  with  accretion  luminosity.  Furthermore,  Najita  et  al. 
(2007)  report  variability  in  BP  Tau’s  short  wavelength  NIR  ex¬ 
cess.  Dutrey  et  al.  (2003)  suggest  that  BP  Tau  may  represent 
a  class  of  transition  object  in  the  process  of  clearing  out  its 
disk  based  on  anomalously  low  12CO  to  dust  emission  ratios. 
Dutrey  et  al.’s  imaging  of  CO  emission  raises  the  possibility  that 
the  CO  photospheric  lines  used  to  measure  the  stellar  RV  may 
be  contaminated  by  variable  CO  emission  caused  by  asymme¬ 
tries  in  a  warm,  rotating  disk.  A  disk  or  accretion  component 
to  the  spectrum  modeling  may  therefore  be  necessary  to  fully 
disentangle  the  various  possible  sources  of  k'-band  variability. 
In  future  work,  we  plan  to  combine  extant  photometric  data 
from  the  McDonald  Observatory  0.8  m  with  our  spectra  to  help 
clarify  the  underlying  source  of  RV  variability  in  not  only  BP 
Tau  but  all  of  our  targets. 

Cl  Tau,  DN  Tau,  and  possibly  V836  Tau  also  illustrate  the 
importance  of  multiwavelength  observations  in  the  case  of 
heavily  spotted  stars.  In  these  cases,  the  optical  spot  noise 
completely  masks  what  may  be  smaller  RV  variability  from  an 
orbiting  planet.  Cl  Tau  is  particularly  interesting  in  this  regard: 
the  strong  correlation  between  optical  RVs  and  bisector  spans 
( r  —  0.7)  leads  to  the  conclusion  that  the  optical  modulation 
is  entirely  activity  induced.  Nevertheless,  we  see  strong  NIR 
variability  at  a  different  period  than  that  seen  in  the  optical  data. 
Furthermore,  subtracting  the  NIR  variability  from  the  optical 
RVs  reveals  a  significantly  cleaner  modulation.  In  this  case,  we 
have  a  target  where  the  optical  RVs  are  almost  entirely  spot 
dominated  while  the  NIR  RV  variations  have  a  different  origin. 
Flere,  our  approach  of  filtering  targets  by  their  optical  activity 
fails.  A  far  more  preferable  approach  is  to  start  with  NIR  RV 
monitoring.  Unfortunately,  we  are  limited  by  practicality.  A 
lack  of  community-accessible,  efficient,  high-resolution  NIR 
spectrographs  dictates  our  emphasis  on  optical  RVs  for  target 
selection. 

Additionally,  we  tested  the  significance  of  the  periods  for  all 
targets  with  Monte  Carlo  simulations.  In  general,  we  found  that 


the  targets  with  the  greatest  number  of  observations  (e.g.,  BP 
Tau,  Hubble  I  4)  have  reasonably  low  (<6%)  FAPs.  For  most 
other  targets,  the  low  number  of  observations  combined  with 
a  high  degree  of  variability  make  the  periods  less  certain.  To 
more  firmly  establish  the  true  periods,  and  eventually  determine 
orbital  parameters  for  likely  companions,  many  more  high- 
cadence  NIR  observations  are  needed. 

6.  CONCLUSIONS 

Consistent  monitoring  of  RV  standards  and  known  exoplanet 
hosts  has  demonstrated  the  long-term  viability  of  our  techniques 
at  a  time  when  NIR  RV  observations  are  coming  into  their  own. 
With  a  stability  of  66  m  s'1  over  four  years,  over  20%  of  the  ex¬ 
tant  exoplanet  population  is  accessible  to  our  observing  method¬ 
ology.  Furthermore,  this  success  has  been  demonstrated  on  a  3  m 
class  telescope  using  a  20  year  old  Cassegrain-mounted  spectro¬ 
graph.  While  parallel  efforts  on  large  telescopes  (Keck,  VLT)  are 
highly  encouraging,  success  on  a  small  telescope,  where  there 
is  greater  community  access,  is  arguably  of  even  greater  value. 
The  large  blocks  of  consecutive  nights  required  for  a  project 
like  this  are  very  difficult  to  acquire  on  8-10  m  class  telescopes. 
This  work  supports  the  findings  of  the  recent  ReSTAR  (Renew¬ 
ing  Small  Telescopes  for  Astronomical  Research)  report:  small 
(<6  m)  telescopes  continue  to  produce  innovative  science  that 
not  only  complements  efforts  on  8-10  m  class  facilities  but  also 
explores  niches  for  which  large  telescopes  are  not  well  suited. 

While  we  have  demonstrated  the  capabilities  of  our  tech¬ 
niques,  we  have  yet  to  confirm  the  presence  of  any  planets  in 
our  young  star  sample.  Mayor  et  al.  (2011)  report  that  65%  of 
Sun-like  stars  harbor  a  planet  of  any  mass  with  P  <  100  days. 
Of  the  known  RV-detected  exoplanets  with  P  <  100  days  (such 
that  their  orbits  are  well  sampled  in  our  survey),  4%  have  semi¬ 
amplitudes  greater  than  the  200  m  s-1  needed  for  a  robust  3a 
detection  with  CSHELL.11  Since  the  T  Tauri  stars  of  today 
will  evolve  to  become  the  solar  analogues  of  tomorrow,  we  ex¬ 
pect  3-4  planet  detections  out  of  the  entire  143  star  sample 
(143  x  0.65  x  0.04)  if  the  PMS  planet  population  mirrors  the 
main-sequence  one.12  Currently,  only  Cl  Tau  is  a  strong  con¬ 
tender;  the  periodicity  seen  in  the  other  two  candidates,  DN  Tau 
and  V836  Tau,  is  so  far  only  suggestive. 

The  lack  of  planet  detections  to  date,  while  likely  impacted 
by  small  number  statistics,  may  indicate  a  real  difference 
between  the  main  sequence  and  T  Tauri  planet  populations. 
In  situ  formation  of  hot  Jupiters  is  unlikely  given  both  the 
unreasonable  requirements  for  delivery  of  solids  to  the  inner 
disk  and  accretion  efficiency.  Therefore,  short-period  massive 
planets  most  likely  arrive  at  their  final  orbits  via  migration  (e.g., 
Kley  &  Nelson  2012).  If  hot  Jupiters  are  not  present  around 
T  Tauri  stars,  this  can  place  a  lower  limit  of  1-3  Myr  on  the 
migration  timescale.  In  this  scenario,  the  Jupiters  have  formed, 
most  likely  beyond  the  snowline,  but  have  not  yet  migrated 
into  short-period  orbits.  If  hot  Jupiters  arrive  via  a  migration 
mechanism  with  a  timescale  greater  than  several  million  years, 
a  survey  of  young  star  associations  over  a  range  of  ages  would 
be  beneficial.  By  identifying  the  age  at  which  hot  Jupiters  are 
detected,  such  a  survey  could  place  observational  limits  on  the 
migration  timescale. 


11  http://exoplanets.org 

12  These  numbers  appear  sensitive  to  selection  criteria.  Broadening  our  criteria 
to  planets  with  any  RV  data  (including  RV  follow-up  of  transit  detections)  as 
opposed  to  limiting  the  sample  to  those  systems  discovered  by  RV  surveys 
increases  the  expected  number  of  detections  to  eight  to  nine  planets. 
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After  seven  years  of  observations,  we  have  also  yet  to  discover 
new  short-period  stellar  or  brown  dwarf  companions.  Similar 
findings  have  been  uncovered  for  brown  dwarf  companions 
to  main-sequence  stars.  While  ~16%  of  solar  analogues  have 
close  (P  <  5  years)  companions,  less  than  1%  of  those  are 
brown  dwarfs;  the  remainder  are  stellar  (~11%)  or  giant  planet 
(~5%)  companions  (the  “brown  dwarf  desert”;  e.g.,  Grether 
&  Lineweaver  2006).  Even  before  completion  of  our  survey, 
this  null  result  on  short-period  brown  dwarf  detections  to  date 
suggests  that  the  desert  is  not  evolutionary  in  nature  (e.g., 
Armitage  &  Bonnell  2002)  but  rather  intrinsic  to  the  formation 
of  brown  dwarfs. 

T  Tauri  stars  reside  in  complex  environments  and,  as  such, 
require  careful  observations  and  analyses  to  isolate  the  various 
origins  of  their  observed  variability.  We  may  find  that  between 
star  spots,  changing  effective  temperatures,  shocks,  and  disk 
obscuration,  these  young  stars  are  too  poorly  suited  to  RV 
analysis.  All  other  approaches,  however,  have  just  as  many,  if 
not  more,  drawbacks.  Direct  imaging  is  the  only  technique  that 
is  beginning  to  successfully  sample  the  young  star  population, 
but  it  is  limited  to  very  massive  planets  on  wide  orbits.  Only  an 
RV  survey  can  adequately  sample  the  inner  AUs — possibly  the 
terrestrial  planet  forming  region.  The  potential  reward  therefore 
justifies  the  intensive  effort  required.  While  the  discovery  of  any 
one  planet  is  unlikely  to  discriminate  between  diverse  planet 
formation  models,  the  confirmation  of  a  planet  around  any  of 
our  targets  would  be  paradigm  shifting.  It  would  represent  the 
first  solid  data  point  of  a  hot  Jupiter  from  the  epoch  of  active 
planet  formation  and  provide  one  of  the  first  observational  limits 
on  the  formation  timescale. 
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